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SPECTROSCOPIC CONFIRMATION OF TWO LYMAN BREAK GALAXIES AT REDSHIFT BEYOND 7 
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ABSTRACT 

We report the spectroscopic confirmation of two Lyman break galaxies at redshift > 7. The galaxies 
were observed as part of an ultra-deep spectroscopic campaign with FORS2 at the ESO /VLT for the 
confirmation of z ~ 7 "z-band dropout" candidates selected from our VLT/Hawk-I imaging survey. 
Both galaxies show a prominent emission line at 9735A and 9858A respectively: the lines have fluxes 
of ~ 1.6 — 1.2 x 10 -17 ergs^ 1 cm~ 2 and exhibit a sharp decline on the blue side and a tail on the 
red side. The asymmetry is quantitatively comparable to the observed asymmetry in z ~ 6 Lya 
lines, where absorption by neutral hydrogen in the ICM truncates the blue side of the emission line 
profile. We carefully evaluate the possibility that the galaxies are instead at lower redshift and we are 
observing either [O n] , [O in] or Ha emission: however from the spectroscopic and the photometric 
data we conclude that there are no other plausible identifications, except for Lya at redshift > 7, 
implying that these are two of the most robust redshift determination for galaxies in the reionization 
epoch. Based on their redshifts and broad-band photometry, we derive limits on the star formation 
rate and on the ultraviolet spectral slopes of the two galaxies. We argue that these two galaxies alone 
are unlikely to have ionized the ICM in their surroundings. 

Subject headings: galaxies: distances and redshifts - galaxies: high-redshift - galaxies: formation 



1. INTRODUCTION 

The recent progress in the exploration of the very 
early Universe has been quite remarkable with several 
groups beginning to assemble large samples of candidate 
high-redshift galaxies (z > 6.5). These include narrow- 
band select ed Ly a emitte rs (L AEs) (|Ota et alj 120101 : 
lOuchi et al.ll2010l : lHu et al.ll2010| ). and color-selected Ly- 
man break galaxies (LBGs) from wide and deep sur- 
yeys carried out in the near-IR, mainly with WFC3 
(iBouwens et al.l l2010at iMcLure et al.1 120101) and H awk- 
I (|Castellano et al.ll2010at iCastellano et alj|2010bl ClOa 
and ClOb in the following). 

Despite the large number of candidates, spectroscopic 
confirmation of a sizable sample of high redshift galaxies 
is still lacking. A systematic spectroscopic follow up is of 
paramount importance for three different reasons. First, 
to evaluate the accuracy and reliability of the Lyman- 
break selection technique that, albeit quite successful at 
z ~ 6, has never been validated at these extreme red- 
shifts. In addition, the intensity of the Lya can pro- 
vide clues on the rest frame properties of the observed 
galaxies: knowledge of the exact redshift and the contri- 
bution of Lya and other emission lines to the observed 
broad-band magnitudes are needed to derive the rest- 
frame properties of the galaxies from the observed con- 
tinuum (SFR, dust content, mctallicity). Finally, and 
more intriguing, the very visibility of the Lya and the 
distribution of its equivalent width in high-z galaxies can 
provide useful constrains on the ionization state of the 
intergalactic medium (IGM) at epochs less than 800 Myr 
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after the Big-Bang (e.g. Fontana et al. 2010, F10 in the 
following; Ouchi et al. 2010). 

To date, only a few galaxies have been confirmed at 
z ~ 7 or beyond. Followin g the early ident ification of 
the z = 6.96 Lya emitter by live et all ((2006). a gamma 
ray burst host was reported at z ~ 8.1 by both Sal- 
vaterra et al. (2009) and Tanvir et al. (2009), but with 
considerable uncertainties on the exact redshift, due to 
the absence of an emission line in the spectra. Very re- 
cently the det e ction of a galaxy at z=8.6 was claimed by 
iLehnert et alj (|2010l ). although the Lya emission line in 
the spectrum has very low S/N and a high uncertainty 
in the flux calibration. 

In this context we have started a systematic campaign 
of spectroscopic follow-up of z ~ 7 "z-band dropout" 
candidates, selected from our imaging survey obtained 
with VLT/Hawk-I (ClOa and ClOb). In F10 we pre- 
sented the results on the sample selected in the GOODS- 
S field. Out of seven candidates observed, we tentatively 
detected only one weak Lya emission line at z = 6.97. 
As discussed in that paper (see also Stark et al 2010), 
this very low fraction of confirmations is at odds with 
what is expected by extrapolating the z — 5 — 6 surveys 
which detect a much larger fraction of Lya emitters. 

Our survey has continued over the two other fields de- 
scribed in ClOb. In this paper we present the only two 
objects with a clear Lya line at z > 7 found in the en- 
tire survey, which are both detected in the BDF4 field 
(Lehnert & Bremer 2003). The results on the final sam- 
ple with a full discussion of their implication is deferred 
to a forthcoming paper (Pentericci et al. in preparation) . 

All magnitudes are in the AB system, and we adopt 
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Ho = 70 km/s/Mpc, M = 0.3 and Q A = 0.7. 

2. TARGET SELECTION AND SPECTROSCOPIC 
OBSERVATIONS 

The targets were selected according to the criteria de- 
scribed extensively in ClOb: besides the three z-band 
dropout candidates listed in Table 3 of that paper, other 
slits were filled with less secure z — 7 candidates as well 
as i-band dropouts. 

Observations were taken in service mode with the 
FORS2 spectrograph on the ESO Very Large Telescope, 
during July- August 2010. We used the 600Z holographic 
grating, that provides the highest sensitivity in the range 
8000 - 10000 A with a spectral resolution R ~ 1390 and 
a sampling of 1.6 A per pixel for a l"xl2" slit. The data 
presented here come from the co-addition of 86 spectra 
of 665 seconds of integration each, on a single mask, for a 
total of 15.9 hr, with median seeing around 0.8". Series 
of spectra were taken at two different positions, offset 
by 4" (16 pixels) in the direction perpendicular to the 
dispersion. 

Standard flat-fielding, bias subtr action and wavelength 
calibration have been applied as in IVanzella et all (2009) 
and F10. The sky background has been subtracted be- 
tween consecutive exposures, exploiting the fact that the 
target spectrum is offset due to dithering. Before com- 
bining frames, particular care has been devoted to the 
possible offset along the wavelength direction, by mea- 
suring the centroids of the sky lines in the wavelength 
interval 9400-9900A. We have also carried out the sky 
subtraction by fitting a polynomial function to the back- 
ground. The two approaches provide consistent results. 

Finally, spectra were flux-calibrated using the obser- 
vations of spectrophotometry standards. Slit losses are 
small, given the extremely compact size of the targets 
and have been neglected in the subsequent discussion. 

3. RESULTS 

3.1. Redshift determination 

We detect a prominent emission line in the spectra of 
two galaxies, candidates BDF - 521 and BDF - 3299 
(ClOb), at wavelengths of 9735A and 9858A respectively. 
In Figures Q] and [5] we present the sky-subtracted ex- 
tracted 1-dimensional spectra and 2-dimensional spec- 
tra for both objects. No other lines are detected in 
the rest of the spectra and no continuum is detected 
for either object. The total line fluxes are 1.6±0.16 x 
10~ 17 ergs -1 cm~ 2 and 1.2±0.14 x 10 -17 ergs 1- cmT 2 re- 
spectively. 

Both lines show a clear asymmetric profile with a sharp 
decline on the blue side and a prominent tail on the red 
side. This asymmetric profile, which we attribute to ab- 
sorption by neutral hydrogen, is the best and in many 
cases unique diagnostic of high-z Lya emission. Indeed 
most LAEs and LBGs at high redshift are too faint to 
detect the break in the continuum caused by IGM at- 
tenuation, with the ex ception of few very bright objects 
(|Kodaira et al. II2003D . 

We have investigated the possibility that the lines 
are instead due to other features, such as Ha, H/3, 
[Om]A5007 or the doublet [On]AA3726 - 3729 in the 
spectra of lower redshift objects. If the detected emis- 
sion lines were H/3 or [O in]A5007 at lower redshift (z^l.O 



and ~0.95), then both emission lines should be seen in 
the observed wavelength range. Actually the other com- 
ponent of the [O in] doublet (at A4959) should also be 
observed, although with lower S/N. No other lines are 
detected in the spectra. We also checked that the posi- 
tions of the expected lines were not coincident with any 
bright skylines. 

In case of Ha emission from z ~ 0.50 galaxies, the rest- 
frame EW of the line would be exceeding 300A, a value 
that is very rarely observed in star forming galaxies at 
low redshift (|Salzer et al. I [2005] ). Furthermore, we note 
that these objects would have a relatively bright con- 
tinuum also at wavelengths below l/i, and should there- 
fore be observed in the deep R, I, z-bands, where in- 
stead we set very stringent upper limits for non-detection 
(AB = 28 — 29 depending on the band, see Table 1 of 
ClOb). 

Finally, in the case of [O n] emitters, the resolution of 
our spectra (R=1390) would be enough to distinguish the 
two components of the doublet, which at z = 1.6 are sep- 
arated by 8A. We searched for examples in our masks, 
and we do detect the [O n] doublet in two galaxies at 
z~1.6, one at z^l.5 and two at z<1.3. In particular the 
two [On] at z^l.6 fall in the same wavelength region as 
the lines discussed here and have a lower flux compared 
to our candidates. In all cases we can clearly distinguish 
the two components: examples at redshift 1.3, 1.5 and 
1.6 are shown in the inset of Figure [2] (top right). Fur- 
thermore, we note that typically the [O n] doublet even 
when unresolved produces an opposite asymmetry with 
respect to the Lya shape, because the A3726A compo- 
nent is often weaker than the A3729A (see also Rhoads 
et al. 2003). 

Since the correct identification of the line is a criti- 
cal issue when no other spectral features are present, as 
mentioned above, the asymmetry can be used to distin- 
guish high-z Lya emission from foreground [O n], [O III], 
or Ha emitters (see Stern et al. 2000). To test our 
ability to recognize asymmetric lines with the S/N of 
our spectra, we run Monte Carlo simulations by insert- 
ing two-dimensional lines of different shapes (Gaussian, 
truncated Gaussian with variable width, and [Oil] dou- 
blet with variable ratio of the two components) in the 
original science frames at positions corresponding to the 
observed wavelengths, normalized to the flux of our lines. 
The frames have been matched in seeing and spectral res- 
olution and reduced as the real ones. This was repeated 
for all slits where we have no signal detected. We find 
that an asymmetric line is always recognized as such if 
it has a width > 150/cto/s, and that the [Oil] doublet in 
the 1-D spectra is always resolved into two components. 

We also calculated the weighted skewness S w intro- 
duced by Shimasaku et al. (2006), to quantify the asym- 
metry of the two lines. Lya emission produces always 
a positive S w : they set a conservative value of S w > 3 
to distinguish LAEs from foreground emitters, although 
there could be also LAEs with S w < 3. For our two 
galaxies we find S w > 5 — 8 (depending on the wave- 
length's range assumed for the measurement), very sim- 
ilar to the average S w of z = 6.6 LAEs (S^ =e - 6 ) = 
7.31 ± 1.5lA . 

We conclude that both lines can be safely identified 
with Lya emission: this implies redshifts of 7.008±0.002 
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and 7.109±0.002 for BDF-521 and BDF-3299, respec- 
tively. 

3.2. Continuum fluxes and spectral slopes 

In the following we determine the correct UV contin- 
uum magnitudes and the Lya equivalent widths for each 
object. 

• BDF 521: Figure [T] (top panel) shows the rel- 
ative position of the Lya line and the two filters Y 
and z. The IGM opacity affects only 2.7% of the Y- 
band filter (<0.05 mag for a flat F v spectrum), therefore 
we neglect the IGM attenuation. The observed magni- 
tude is Y = 25.86 ± 0.11, corrected to Y cont = 26.31 
if the Lya contribution is subtracted. We note that 
for this source there is also a 2-sigma detection at mag 
28.00±0.60 in the z-band which is consistent with the 
Lya falling at the edge of the z-band filter. Assuming 
Y con t — 26.31 the rest-frame Lya equivalent width is 
64^g°A (where the errors come from the line flux uncer- 
tainty). 

• BDF-3299: The observed magnitude is Y = 26.15± 
0.14, corrected to Y cont = 26.58 if Lya is subtracted. 
However in this case the IGM attenuation affects 15% of 
the V-band filter, absorbing completely the continuum 
blue- ward of the Lya line (see Figured]). Assuming con- 
servatively a flat spectrum for the continuum, we derive 
Y con t = 26.40 and a rest-frame Lya equivalent width of 

sot^A. 

Adopting these Lya EWs we derived the slopes of the 
UV continuum, fitting to the broad-band photometry a 
set of simple models with flux scaling as F\ oc , which 
include also the Lya emission and IGM absorption. We 
derive the best fit models and the relevant error analysis 
with a standard \ 2 minimization. The results are shown 
in Figure 3, where we plot the best-fitting solutions along 
with the error analysis (see insets). The best-fit slopes 
are remarkably steep, with f3 = —3.7 and —2.7 for BDF- 
521 and BDF-3299, respectively, similar to those recently 
claimed for z-band dropouts on the basis of new HST 
WFC3 data (j3 < -2.5, Bouwens et al. 2010b; Oesch et 
al. 2010; Bunker et al. 2010; Finkelstein et al. 2010). 
However, much shallower solutions with /3 < —1.7 and < 
—0.9 are still allowed at 1-sigma c.l. (Ax 2 = 1). Clearly, 
these results are mostly constrained by the non-detection 
in the J and in the K-band (see Table 1). 

3.3. Implication for the stellar population 

We have first explored whether these lines can be as- 
cribed to AGN. From the non-detection of NV in the 
stacked spectrum we derive a ler lower limit on the 
Lya/NV > 8. This value is higher than the average ra- 
tio for QSOs but not unusual for weak narrow-line AGN 
found amongst LBGs (e.g. Steidel et al. 2002), so we 
cannot exclude this possibility. However we rely on the 
fact that AGN are much rarer than galaxies, and in the 
following we will assume that the objects are star forming 
galaxies. 

At the redshifts derived above the implied luminosi- 
ties for the Lya lines are 9.3 and 7.3 x 10 42 ergs _1 . We 
estimate the star formation rate (SFR) from the Lya lu- 
minosity using the Kennicutt's relation (Kennicutt 1998) 
with the case B recombination theory as SFR = 9.1 x 
L(Lya), with L(Lya) in units of 10 43 ergs -1 . We obtain 



values around 7-9 M yr _1 (see Table 1): they represent 
lower limits since they are not corrected for absorption 
effects which depend on various parameters, including 
the neutral fraction of the IGM and the kinematic status 
of neutral hydrogen within the galaxies (e.g. Ahn 2004). 

From the continuum luminosity we obtain 
L(UV 1275 ) = 8.1 x 10 28 and L(UV 12e o) = 
7.6 x lO^ergs^ 1 Hz^ 1 . We convert luminosities 
into SFR using the Kennicutt relation for UV contin- 
uum, SFR{UV) = 1.4 x 1Q- 28 L(UV). We determine 
the flux at 1500A assuming the slope beta derived 
before. We obtain SFRs of 8.9 Mgyr" 1 and 9.4 M^yr' 1 
respectively, similar to those determined from the 
Lya line. These values are obtained assuming no 
dust and the agreement of SFR(UV) and SFR(iya) is 
consistent with this scenario. Using instead a global 
fit to the overall photometry with Bruzual & Chariot 
templates and a Calzetti (2000) attenuation law, we find 
best fit values for SFR that are in excellent agreement 
with the above determinations, and E(B-V)=0. The 
standard error analysis on these fits provides the limits 
one can place on the maximum SFR. Given the poor 
constraints on the J and K-bands, we find that for all 
models acceptable at 68% confidence level E{B — V) 
is lower than 0.25-0.35, while the SFRs upper limits 
are 160M Q ?/r -:L and 400M Q ?/r _:L , respectively, for the 
two objects, depending on the combination of age and 
metallicity. Note that the adoption of a Gallerani et 
al. (2010) attenuation curve would provide a SFR limit 
about 2 times lower. 

4. DISCUSSION 

With the spectroscopic identification of these two z > 7 
galaxies we have moved closer to the epoch of reionization 
and may possibly be observing "re-ionizers" at work at 
the end of this process. While it is clearly impossible to 
draw any statistical result out of two objects only, we 
discuss here a few implications. 

Recently, several authors have argued that the fraction 
of galaxies with large Lya EW increases in the redshift 
range 4 < z < 6 (e.g. Stark et al. 2010; Douglas et 
al. 2010; V09). The two emitters reported here with 
Lya EW>50A seem to continue this trend. Indeed the 
fraction of z-band dropouts spectroscopically confirmed 
is much higher for this field (two out of three candidates 
observed), compared to the GOODS-south field (one out 
of seven candidates, F10). 

Therefore the main results of F10, namely that the 
fraction of galaxies with Lya EW exceeding 50A seems to 
show a sharp decline from z ~ 6 to z ~ 7, thus reversing 
the increasing trend from z ~ 4 to z ~ 6, appears in 
conflict with the present observations. 

However we anticipate that the result in F10 is con- 
firmed by the analysis of the final sample that includes a 
third field (Pentericci et al. in preparation) and a much 
larger statistics. The present spectra indicate that there 
are considerable field to field variations and therefore 
solid conclusions can be drawn only from large samples. 

Turning to the nature of these objects, one can deter- 
mine whether they are able to ionize the IGM around 
them. Many models indeed predict that only intrinsi- 
cally UV bright galaxies have enough photons to build 
reionized "bubbles" around them, on reasonably short 
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timescales (|Daval et al.l l2008f ). Following Loeb et al. 
(2005) and assuming SFR of lOMgyr" 1 and age lOOMyr 
we derive a maximum radius of the H n region of 
Rmax = 0.73(/ esc ) 1/3 Mpc physical with f esc the frac- 
tion of escaping ionizing photons. We are assuming iso- 
lated galaxies surrounded by a medium that is mostly 
neutral when they started to form stars (z ~ 8). In 
order for the Lya to be transmitted the optical depth 
Tdamp to Lya absorption at the galaxy redshift has to be 
less than one (Wyithe & Loeb 2005): this corresponds 
to a R = 1.1 Mpc. Therefore the radius we have derived 
would not be large enough even if a maximum / esc = 1 
is considered. 

As already mentioned the intrinsic SFRs could be con- 
siderably higher: for the maximum SFR consistent with 
our photometry, the galaxies could ionize a large enough 
region if they had f esc > 0.1. Such values of f esc are 
not commonly observed in the z ~ 3 Universe with few 
exceptions, such as the compact LBG at z = 3.8 recently 
reported by Vanzella et al. (2010) for which the Lyman 
continuum (A < 912A) has been detected directly and 
whose non-ionizing features reveals j3 = —2.1, very weak 
interstellar absorption lines and absence of Lya emission. 

We stress that with the current depth of the near in- 
frared observations, it is not possible to produce better 
constraints on the UV slopes, dust content and SFR. This 
issue could easily be solved with a modest investment of 
HST time, by obtaining deep Y, J and H-band WFC3 
observations. 

If the steep slopes and the dust-free scenarios are con- 
firmed, then the galaxies do not have enough photons to 
ionize the surrounding H II regions. In such a case, the 
visibility of the Lya line in these two objects implies the 
existence of additional ionizing sources, either galaxies 
and/or objects of different nature. It is intriguing to sur- 
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mise that the field-to-field variations in Lya visibility we 
observe indicate that the reionization process was very 
inhomogeneous at these epochs. A possible reason is a 
clustering effect, particularly efficient since the first ob- 
jects formed in highly biased regions (e.g. Furlanetto et 
al. 2006). 

We point out that the two confirmed LBGs in the BDF 
field are at a physical distance from each other that is 
of the same order as the clustering length of LAEs at 
z = 6.6 r = 2 - 5hiQ Mpc (Ouchi et al. 2010) (they 
are separated by 6' in the plane of the sky corresponding 
to 1.9 Mpc and by 0.101 in redshift space corresponding 
to a proper distance of 4.4 Mpc). However there is no 
evidence of significant enhancement in the density of z- 
band dropouts in this field. 

Another possible mechanism suggested by Wyithe & 
Loeb (2005) is related to a previous QSO activity, i.e., 
the Hu regions generated by quasars remain as a fossil 
after the quasar activity ends, since the recombination 
time is longer than the Hubble time at the mean IGM 
density for z < 8. 

These scenarios and their combinations can be tested 
only when a sizable sample of high-z sources are con- 
firmed and their clustering and physical properties are 
determined. 
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Fig. 1. — 1-dimensional spectra of candidates BDF-521 and BDF-3299. On the left panels spectra are shown with superimposed the 
z-band and Y— band filters. In the bottom the 1-dimcnsional (flux calibrated) spectrum of the sky is shown (the position of the Lyo lines 
are marked with triangles). On the right side the zoomed Lya lines are shown with the position of the peak marked with a vertical dotted 
line. The two lines are significant at around 20 sigma level. In all panels units in the Y-axis are 10 -18 erg s~ 1 cm~ 2 A - 1 . The red dashed 
lines superimposed to the zoomed Lyct (on the right) show the spectrum of the sky in arbitrary units. 
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Fig. 2. — 2-dimensional sky subtracted spectra of BDF-521 (top) and BDF-3299 (bottom). In both cases the emission line is clearly 
detected. The "negative" features above and below the emission (typical of the sky subtraction technique) are also evident. In the middle 
the sky spectrum is shown as a reference for sky line positions. The resolved [On] doublets at redshift 1.3, 1.5 and 1.6 are shown in the 
upper right corner, from bottom to top, respectively. 
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Fig. 3. — Observed photometry and best-fitting spectral model for BDF-521 (upper panel) and BDF-3299 (lower panel). Black dots and 
upper limits correspond to observed magnitudes. The solid line is a best-fitting model. Red dots are the predicted magnitudes from the 
model. The blue thin line is the dust-extincted Bruzual & Chariot model which provides the maximum SFR with an acceptable x 2 . In the 
inset the x 2 as a function of fi is shown (the red horizontal line marks the Ax 2 = 1). 
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TABLE 1 

Spectroscopic and photometric properties of the two redshift seven galaxies 



ID 


RA, DEC 
J2000 


Redshift 


i(Lya) 


SFR(Lya ) 

Mqvt- 1 


EWrest 

A 


FWHM" 


S w 

A 


z 


Y 


J 

la 


K 
Ict 


BDF-521 
BDF-3299 


336.9444, -35.1188 
337.0511, -35.1665 


7.008±0.002 
7.109±0.002 


1.62±0.16 
1.21±0.14 


8.5 
6.6 


64 
50 


240 
200 


8.2±l 
> 5.9 b 


28.00 
>28.55 


25.86 
26.15 


>26.5 
>26.5 


>26.0 
>26.0 



Note. — f(Lya) in units of 10 17 ergs 1 cm 2 A 1 . (a) Measured directly from the line and corrected for instrumental broadening, (b) This 
considered a lower limit since the red side of the emission line is truncated by the presence of a bright skyline (see Figure 1). 



